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A comparison is given of the various recently pubUshed extractions of the Sivers 
functions from the HERMES and COMPASS data on single-transverse spin asym- 
metries in semi-inclusive deeply inelastic scattering. 

1. Introduction 

Single-spin asymmetries (SSA) in semi-inclusive deeply inelastic scattering 
(SIDIS) off transversely polarized nucleon targets have been under intense 
experimental investigation over the past few years. Substantial asymme- 
tries have been reported in some cases, in particular, with best statistics, 
by the HERMES collaboration for scattering off a proton target. 

The importance of SSA lies in the fact that they provide new insights 
into QCD and nucleon structure. ^"^'^ For instance, the asymmetry in SIDIS 
may contain an angular dependence of the form sin(^ — t/fg), where </) and 
(j)s denote respectively the azimuthal angles of the produced hadron and 
the target polarization vector with respect to the axis defined by the hard 
virtual photon. This angular dependence arises from the so-called Sivers 
effect^ tightly related to notions of an intrinsic asymmetry in the parton 
transverse momentum distribution and angular momenta. Factorization 
theorems"*^^"^^ proven to leading power in the photon virtuality Q provide 
the basis for a QCD description of the process, and allow to extract the 
Sivers function from SIDIS data^""' and to use it for predictions for the SSA 
in the Drell-Yan (DY) process, hopefully to be explored experimentally at 
RHIC, COMPASS and the GSI. Comparisons of SIDIS and the DY process 
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will be particularly important for testing our understanding of the underly- 
ing physics, since it has been predicted^'^'^'' that the Sivers functions appear 
with opposite signs in these two processes. The approach just outlined has 
been followed recently in Refs. [18-23]. In this note we compare the results 
of these papers for the extracted Sivers functions 



In the extractions of the Sivers functions from SIDIS several simplifying 
approximations were common between the groups, namely the neglect of 

the so-called "soft factor" -'^^'-'^^ and the Sivers antiquark functions. Different 
approaches were, however, followed in Refs. [19-23] concerning the treat- 
ment of the dependence of the distributions on transverse parton momenta. 

The Sivers SSA is obtained^'^^ by weighting the events entering the spin 
asymmetry with sin((/)— 0s). When analyzed in this way, however, specific 
models for the dependence on parton transverse momenta need to be made 
in the theoretical expression. By assuming that the transverse momentum 
dependence of the Sivers function is of the form f^{x, p^) = ,f^{x) G(p|,) 
and/or similarly for other distribution or fragmentation functions, the 
Sivers SSA as defined at HERMES^ can be written generically as 

The factor (—2) is due to conventions^^ and Fg^^{x) is some functional 
depending on and the model used for parton transverse momenta. 

Notice that by including in addition a factor of Ph± /Mm into the weight 
in (2) the resulting SSA can be interpreted model-independently in terms 
of the transverse moment of the Sivers function^ ^ 

ft^'^\x) ^ J d'pT ^ AV(x, pI) = - / d^PT ^ A^/,/,r {X, p|) . 

(3) 

Such weighted SSA were argued to be less sensitive to Sudakov suppression 
which can be important for predictions involving the Sivers function. Pre- 
liminary HERMES data for such SSA are available^ and were studied in 
Ref. [18], where a first fit for the transverse moment of the Sivers function 
(3) was obtained. The result of [18] is in good agreement with the stud- 
ies of SSA analyzed without a power of Ph± in the weight^"** reported in 
Refs. [19-23]. The next Sections review and compare the fit results for the 
Sivers functions extracted in the different approaches in Refs. [19-23]. 
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2. The approach of Refs. [19,20] 

In Ref. [19] the azimuthal angular dependence (Cahn effect) of the SIDIS 
unpolarizod cross section was used to extract the widths of the Gaussian 
PT-dependent parton distribution (pdf) and fragmentation (ff) functions 
respectively as = 0.25 (GeV/c)^ and (K^) = 0.2 {GeV/cf. A first es- 
timate of the Sivcrs functions was then obtained by fitting the data on 
^sv^4>-<Ps) observed by HERMES collaboration.^'^ In Ref. [20] a novel 
fit on the new HERMES data^ together with data from the COMPASS 
collaboration^ was performed. In both fits the full exact kinematics was 
always adopted. The Sivcrs function (u, d quarks) was parameterized as: 

A^/g/pT {x, p2 ) = 2 N,{x) f,/p{x) g{p^r) Kp^t) , (4) 

-i-h \{aq+bq) p-Pt/{Pt) 

N,{x) = N.x^^il - xP l^^±Mj^ , g{p%) = . (5) 

Two options for the /i(Pr) function were considered, namely: 

the latter allowing, at leading order in pr/Q, to give for Fg;^ in Eq. (2): 



-Psiv(a;) = V / ^ = /it (2^) with {p\) 



(7) 

In the fits, ,fq/p{x) was taken from the LO MRSTOl set^^, whereas Kretzer's 
set^^ for the LO ff was used. The 7 parameters were then extracted as^°: 

7V„ = 0.32 ±0.11 a„ = 0.29 ±0.35 6„ = 0.53 ±3.58 
A^d = -1.0 ±0.12 ad = 1.16 ±0.47 6,; = 3.77 ±2.59 (8) 
M'2 = 0.55 ± 0.38 (Mo^ = 0.32 ± 0.25) (GcV/c)^ 

with a per degree of freedom (x^^f ) of 1.06. The one-sigma band shown 
in Fig. 1 (Eq. (6b)) takes into account the errors with their correlations. 

These results were then used to give predictions for SSA measurable in 
SIDIS and DY processes for various kinematical configurations. 

These effects were also invoked^'^"'^^'^^ to generate SSA for other pro- 
cesses in hadron-hadron-collisions^^'"^" although the status of factorization 
is less clear in this case. Here we only point out that the SIDIS data are 
sensitive to much smaller x values than the E704 (STAR) ones. 
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3. The approach of Ref. [21] 

In Ref. [21] it was assumed that the final hadron's transverse momentum 
is entirely due to the transverse-momentum dependence in the Sivers func- 
tion. There is then no further assumption on the particular form of this 
dependence; rather it is integrated out in order to compare to the exper- 
imental data. The transverse momenta contributed by the other factors 
in the factorization formula will give some smearing effects which may be 
viewed as "sub-dominant". (However, we emphasize that this will not be 
true toward small z where the transverse momentum in the fragmentation 
functions will become important, likely resulting in a suppression of the 
asymmetry at small z.) The "1/2-moments" of the Sivers functions were 
then introduced in Ref. [21] in the fit to the experimental data: 

/rfV^A"T^(x,p|). (9) 

These appear in an expression of the form (2) for the Sivers asymmetry, 
where 

^liv(^) = ^ • (10) 

In the actual fit to the HERMES data in^^ the functions q"^^"^^ [x) were 
modeled in terms of the unpolarized w-quark distribution as 

-Sux{l-x) , ^^-P. = - , (11) 



u{x) u{x) 

where u{x) was taken from the GRV LO parameterizations for the unpo- 
larized parton distributions.^^ Furthermore, Kretzer's set for the LO frag- 
mentation functions^^ was used. The fit to the new preliminary HERMES 
data gave 

5„ = -0.81 ± 0.07, S'd = 1.86 ±0.28 , (12) 

with Xdof ~ A fit to the old published HERMES data gave instead Su = 
-0.55 ±0.37 and 5^ = 1.1± 1.6, with a similar size of Xdof ■ The COMPASS 
data were not included in the fit performed in^^, but a comparison of the 
fit with the data was given, showing good agreement. The results of the 
fit to the HERMES data were furthermore used for making predictions for 
the SSAs in the Drell-Yan process and in di-jet and jet-photon correlations 
at RHIC. 
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4. The approach of Refs. [22,23] 

In Rcf. [22, 23] the distributions of transverse parton momenta in 
and Di were assumed to be Gaussian with the respective widths 
(PT)siv and (K^) taken to be flavour- and x- or z-independent. In this 
model the F^^^ defined in (2) is given by the expression in Eq. (7) with 

(pI) replaced by {p'^}sw- 

The values {K^) = 0.16 (GcV/c)^. (p^) = 0.33(GeV/c)2 were 
extracted^^ fj.^^ ^j^g HERMES data'^^ (^p^^j .^^^ 

are similar to those 

discussed in Sec. 2, while {'p^)siv S [0.01; 0.32] (GeV/c)^ remained poorly 
constrained by positivity'^'' still allowing an extraction of the transverse 
moment of the Sivers function (3). 

In order to reduce the number of fit parameters the prediction^^ from 
the limit of a large number of colours Nc was imposed: 

/fi^(a;,p|n) = — /f^(x, py) modulo 1/A^c corrections. (13) 

The best fit^^ (using parameterizations^^'^®) to the published data? is 

xftr^'^ix) " Ax\l - xf -0.17a;0-66(l - xf (14) 

with a Xdof ~ O-"^' ^^'^ ^ ^''^ uncertainty of roughly ±30%. This re- 
sult agrees well with the fit to the preliminary P/ij^-weighted HERMES 
data^, which were analyzed in a (transverse parton momentum) model- 
independent way [18]. The good agreement of the results in Refs. [18, 22] 
is an important cross check for the applicability of the Gauss model to the 
description of SSA in SIDIS. 

For sake of a better comparison to the results by the other groups^^'^°'^^ 
the above fit procedure was applied^^ to the most recent and more precise 
preliminary HERMES data.^ The new fit has a Xdof ^ ^ ^"^^ ^ consistent^^ 
with that quoted in Eq. (14). One has to keep in mind that the large-A^c 
relation (13) is a useful constraint at the present stage, and will have to be 
relaxed when future more precise data will become available. 

Note that for (K^) ^ in (7) one obtains Fg^^ix) yiT^^^^^ix) 
within the Gaussian model. This limit means that the produced hadron ac- 
quires no additional transverse momentum from the fragmentation process, 
i.e. £>?(^,K|,) = £)°(2)5(2)(Kt). In this sense, the approach of Ref. [21] 
discussed in Sec. 3, c.f. Eq. (10), is contained as a limiting case in the 
Gauss ansatz. 
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Figure 1. The first and 1/2-trans verse moments of the Sivers quark distribution func- 
tions, defined in Eqs. (3, 9), as extracted in Refs. [20, 21, 23]. The fits were constrained 
mainly (or solely) by the preliminajry HERMES data* in the indicated x-range. The 
curves indicate the l-a regions of the various pajrameterizations. 



5. Comparison of the results and Conclusions 

It should be stressed that the various fit results, when used within the 
respective approaches, provide equally good descriptions of the HERMES 
and COMPASS data. Here we compare only those analyses^^'^^'^^ in which 
the most recent and more precise preliminary HERMES data^ were used. 

In Fig. la we compare the fits for f^rj;^'''^ from Refs. [20, 23], and in 
Fig. lb the fits for f^r^^^^^'' from Refs. [20, 21]. (A direct comparison of 
[21] and [23] is not possible.) In view of the different models assumed for 
the transverse parton momenta and the varying fit Ansatze, we observe 
a satisfactory qualitative agreement — in the cc-region constrained by the 
HERMES data. However, a closer look reveals differences between the 
results in Fig. 1, which indicate the size of the systematic uncertainties of 
the three Sivers function fits mainly due to the use of different models for 
the parton transverse momenta. These uncertainties were not estimated in 
Refs. [20, 21, 23]. 

We have presented a comparison of three cxtractions^'^'^-'^'^'^ of Sivers 
functions from HERMES and COMPASS data on single-transverse spin 
asymmetries in SIDIS. The three approaches somewhat differ, but they 
describe the data with similar quality. The fits are in good qualitative 
agreement, though there are differences with regard to the size and shape 
of the extracted Sivers functions. These differences reflect the model depen- 
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dence of the fit results which gives rise to a certain theoretical systematic 
uncertainty of the fit results. The latter seems, however, less dominant 
than the statistical uncertainty of the fits at the present stage. 

It is clear that further information from experiment will be vital. For 
now, one cannot really expect to obtain much more than a first qualitative 
picture of the Sivers functions. We also emphasize that it will be crucial 
for the future to experimentally confirm the leading-power nature of the 
observed spin asymmetries. For this, forthcoming COMPASS or JLab data 
for scattering off a proton target and studies of the Q^-dependence of the 
asymmetries will be important. 

The good qualitative agreement between the different approaches ob- 
served here means that the prcdictions^^"^'' for the magnitude of the Sivers 
effect in DY are robust — in the kinematic region constrained by the HER- 
MES data. This solidifies the conclusions^*"^^ that the predicted sign rever- 
sal of the Sivers function between SIDIS and DY, can be tested in running 
or future experiments at RHIC, COMPASS and PAX. 
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